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FORMULAS  OF  TML  SHOH  J’ VV  A V 1'.  ASYMFTO  JIC 
IN  A  PROBLEM  OF  DIFFRAC'IION  1)Y  CONVEX  1K)DIES 

by 

V.  S.  Busl.u  v 


§  ^  •  INTRODUC'J  TON 

1.  i'uriiiulus  are  derived  in  artiile  iur  the  she rtwavt: 

asymptotic  of  Green’s  function  G(x,  x’;  k)  of  tlu;  external  Dirichlet 
problem  for  the  Helmholtz  ecpiation  in  a  plam  .  GveenCs  function 
satisfies  the  equation 

( —  vi  —  A’)  U  ( A-,  a' ;  k)  i( A  A  )  {X,  x'(:D,  k  >  0 )  ^  1  ^ 


and  the  boundary  conditions 


0{x,  a';  k) 


0. 


4  «  L 


j " 

'A* 


<'1/  (X,  a';  k)  j  ,  C..  r, 

- a, A, - X',  A*) -->0, 


(1.2) 


Region  D,  in  which  the  equaticui  is  i  xamincd,  is  outside  a 
finite  closed  convex  contour  L.  It  is  assumed  that  the  radius  of 
curvature  o(s)  of  this  contour,  as  a  fuiu  iiun  -af  the  arc  length  s, 
has  two  continuous  derivatives  <ind  p(s)  :*  >  0. 


Figure  ! 


•  2  « 


The  following  definitions  will  be  used  in  connection  with  the 
geometric  characteristics  of  the  shortwave  asymptotic  (the  asymp¬ 
totic  when  k  f  oo).  Relative  to  point  x',  region  O  is  divided  into 
two  parts:  the  illuminated  part  (light)  R(x')  and  the  shadow  part 
S(x').  The  tangents  H^T  to  the  contour  L  (figure  1), 

whose  extensions  pass  through  the  point  x'  are  the  boundaries  (geo¬ 
metric)  of  these  regions. 

2.  Many  papers  have  been  written  on  the  derivation  of  formulas 
of  the  shortwave  asymptotic  in  which  asymptotic  formulas  are  obtain¬ 
ed  for  a  considerable  part  of  the  characteristic  locations  of  the  points 
X  and  x'  in  the  region  D.  However,  even  in  the  simplest  cases,  there 
is  no  exact  justification  of  these  formulas.  In  an  attempt  to  make 
such  a  justification,  by  using  a  certain  procedure,  it  was  necessary 
to  refine  the  asymptotic  formulas  in  specific  respects.  This  refine¬ 
ment  is  also  described  in  this  article. 

Keep  the  following  system  of  justification  in  mind: 

Let  the  asymptotic  Q(x,  x*;  k)  of  Green's  function  G(x,  x';k),  which 
satisfies  the  conditions  (1.2),  be  known  from  any  non-rigorous  argu¬ 
ments  with  any  x  and  x'(x,  x'  ^  D).  This  residual  is 

Kix,  X';  /:)«.(-v=-A'»)Q(x.  jc';  k)^l{x-x‘)  (1.3) 

and  the  relationship 

Q{x',  x\  A)  =  C/(x',  x\  A)-|-  f  </yO'(x',  y;  k)K(y,  x\  k) 


(1.4) 


ii  conniderud  as  iIr  .*quiituiii  iv.>r  int!  function  G(x’,  x;  k)  (x’  being 
fixed)*  If  the  asymptotic  Q(x,  k)  is  such  that  the  corresponding 
residual  K(x,  x';  k)  generates  the  operator  K  with  a  norm  which 
approaches  zero  as  k  -►  ae  in  some  proper  functiotial  space,  equa¬ 
tion  (1.4)  gives  an  estimate  of  the  error  in  the  asymptotic  Q(x,  x*;  k). 

Familiar  asymptotic  formulas,  which  can  be  found,  e,g, ,  in 
[1]  and  [2-5],  indicate  well  the  overall  structure  of  an  asymptotic 
but  we  have  been  unable  to  construct  that  functional  space  in  which 
the  residual  corresponding  to  these  formulas  generates  an  operator 
with  a  small  norm  (the  asymptotic  Q(x,  x;  k)  should  also  belong  to 
this  space). 

Section  2  of  this  work  gives  a  construction  of  the  asymptotic 
Q(x,  x':  k),  with  the  following  properties: 

a)  Q(x,  x';  k)  satisfies  conditions  (1*2); 

b)  Q(x,  x';  k)  5  Q(x’,  x;  k); 

c)  Q(x,  x'i  k)  is  a  continuous  function  of  the  arguments  x  and 
x',  except  at  the  point  x  r  x',  where  it  has  a  singularity 
characteristic  for  Green's  function; 

d)  Q(x,  x';  k)  has  continuous  second  order  derivatives  with 
respect  to  the  variable  x,  except  at  the  geometric  light- 
shadow  boundary,  where  the  derivative  may  have  discon¬ 
tinuities  of  the  first  kind. 

These  properties  make  it  possible  to  consider  equation  (1,4) 
in  some  space  of  continuous  functions.  A  total  estimate  of  the  norm 
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01  the  operator  K  will  be  treated  in  another  work;  here  we  restrict 
ourselves  only  to  a  partial  study  of  the  properties  of  the  residual 
K(xt  x';  k):  let  us  estimate  its  order  when  k  -f  oo  .  In  section  3 
it  is  shown  that  the  main  terms  of  the  residual  for  the  obtained 
asymptotic  formulas  are  reduced  when  k  -f  oo , 

The  formulas  given  in  section  Z  arc  directly  connected  with  the 
results  in  works  [2-5J.  These  results  are  characterized  by  a  special 
type  of  contour  integrals  by  which  the  asymptotic  is  described.  V.  A. 
Fok  was  the  first  to  introduce  and  investigate  such  integrals. 

$  2.  CONSTRUCTION  OF  THE  ASYMPTOTICS 

1.  Formulas  which  are  known  in  different  cases  (for  special 
contours,  admitting  the  separation  of  variables)  indicate  that  the 
asymptotic  has  the  form  of  series  2^  .  F  (<)>)•  to  each  of  whose  terms 
corresponds  a  specific  phase  The  phases  can  be  described  as 
the  lengths  of  certain  extremum  lines,  connecting  the  points  x,  x' 
and  lying  in  region  D,  If  the  points  x  and  x'  are  located  in  the 
shadow  relative  to  each  other,  they  are  the  shortest  curves,  envelop¬ 
ing  the  contour  L  (figure  1)  and  an  infinite  number  of  smooth  curves, 
which  differ  from  those  shown  by  additional  turns  around  the  contour  L, 
All  these  curves,  the  corresponding  phases  and  components  in  the 
asymptotic,  are  calledenvelopingrays,  phases  and  waves.  When  the 
point  X  approaches  the  geometric  light-shadow  boundary  and  passes 
into  the  illuminated  part,  one  of  the  enveloping  phases  ceases  to  have 


y 


«  . 
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meaning  and  ii  replaced  by  two  phases  (figure  2):  *  |x  •  x'|  and 

which  is  the  length  of  the  broken  line,  reflected  from  the  contour, 
according  to  the  laws  of  geometric  optics  (the  angle  of  reflection  is 
equal  to  the  angle  of  incidence).  The  remaining  enveloping  phases, 
including  those  which  were  obtained  from  the  phases  "having  been  split" 
by  the  additional  terms,  are  now  preserved.  Here,  the  designations 
direct  (for  ^  )  and  reflected  (for  rays  (phases,  waves)  are  used. 
Thus,  there  are  three  kinds  of  waves:  one  direct,  one  reflected,  defined 
for  x6R(x')  and  an  infinite  number  of  enveloping  waves.  For  contrast 
to  and  any  of  the  enveloping  phases  will  be  designated  by 
All  phases  introduced  satisfy  the  equation  of  the  eikonal  (x^  D) 

(2.1) 

The  component  F  (4>)  in  the  asymptotic  has  the  following  struc¬ 
ture: 

r(q.)  *).  (2  2) 

i.  the  Hanktil  function  (  6] .  x';  k)  is  a.eumed  to  be 

0 

a  slowly  varying  function  in  the  following  sense  when  k  f  oo  : 

A  ^  (k<i>)  -  k  Nk<t>)  has  an  order  of  growth  kH^  ^  (k4>).  when 

X  0  ■  v„  0  0 

k  «  the  gradient  V  (Xi  x*;  k)  has  an  order  less  than 

kU^*^^  (x,  x'j  k)  as  k  . 

Let  us  assume 


(A.  jc';  A)  -  -l, 


U.i) 


*  (> 


MO  that  llu  WmV».'  in  the  li^hl  nerves,  with  tlu- 

Mohitiun  of  equation  (1.  1). 


Let  ut  designate  {x»  x';  k)  ■  u'^R^x.  x';  k)  and 

(X,  X';  k)  ■  (x,  X';  k). 

Considering  the  described  structure  of  the  asymptotic  as  the 

(R) 

assumption,  let  us  look  for  the  functions  U'  '  (x,  x';  k)  and 
(SI 

U'  '  (x,  x';  k).  Taking  into  account  that  the  phases,  which  generate 
the  rapidly  variable  factors,  are  different,  it  will  be  required  that 
each  of  the  components  F (^)  asymptotically  satisfy  the  Helmholts 
equation  (for  those  values  of  x  and  x'  when  the  corresponding  com¬ 
ponents  are  defined).  From  those  same  considerations  let  us  require 
that  each  component  satisfy  the  boundary  condition  (1.2)  on  the  contour, 
with  the  exception  of  the  direct  and  reflected  waves  (on  the  contour 
L.  ♦o  *  ♦«).  for  which  we  require  that  their  sum  satisfy  the  boundary 


condition 


-  H  • 


i  ( M 

!)  j.iiK t'.r.  llaiiki  1  luiu  tiDn  H'  '  ii»  uscci  and 

4  0 

t  ijn.i  t  u  1 M  {d  .  I ) , 

.'iiiii  e  tlur  lunction  U(x,  x' ,  k)  is  ciiisunicd  to  change  slowly,  in 
uhlaiuing  t)if  asymptotic  equation  lor  U  (x,  x' ;  k),  we  should  use  the 
asymptotic  form  of  the  Hankel  luiu  tion  whenk4)-^+  oo;  therefore  it 
is  natural  to  preserve  asymptotically  only  the  second  comjionent  in 
the  last  expression  (2.6),  since  both  components  in  the  brackets  con¬ 
tain  the  factors  //*,”' (A-<I))—  o  *)and  O  f(A(I>)  •) 

respectively.  In  this  section,  let  us  agree  not  to  differentiate  between 
the  asymptotic  and  exact  equalities  in  the  notations.  This  difference 
will  be  obvious  each  time  for  the  text.  Asymptotically  we  have 

I- — 0.  (2,7) 


If  only  the  term  in  the  brackets  containing  the  factor  k  remains 

in  equation  (2.7)  and  if  it  required  that  the  solution  of  the  obtained 

(R) 

equation  satisfy  the  boundary  conditions  (2.4),  the  functions  U'  '  (x,  x';k] 
(S) 

and  '  (x,  x';  k)  are  determined  uniquely  (explicit  expressions  are 
in  paragraphs  "2"  and  "3"  of  this  section).  It  is  evident  from  the 

(R)  is) 

explicit  expressions  that  '  and  U'  '  do  not  satisfy  condition  (2.  5) 

of  continuous  merger  at  the  light- shadow  boundary.  In  addition,  the 

(R) 

derivative  of  the  function  U'  '  with  respect  to  the  normal  to  this  bound¬ 


ary  thus  obtained  reverts  to  infinity  at  the  boundary.  Due  to  these 
singularities  no  meaning  can  be  given  to  the  operator  K  in  equation 
(1.4),  which  generates  such  an  asymptotic.  This  indicates  that  some 
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componuntd  ot  thv  .Laplacian  operator  should  be  retained  asymptotically 
in  equation  (^.7),  It  ia  natural  to  expect  that  in  the  vicinity  of  the 
geometric  light- shadow  boundary,  there  should  be  the  second  order 
derivatives  with  respect  to  some  direction  orthogonal  to  the  boundary. 
We  will  arrive  at  such  a  conclusion  by  relying  on  the  anctlogy  with 
constructions  of  the  boundary-layer  type  and  by  taking  into  account 
that  the  derivative  with  respect  to  the  tangent  to  the  boundary  is  already 
contained  in  the  brackets  with  the  factor  k.  From  these  considerations, 
only  the  second  order  derivative  with  respect  to  the  direction  orotho- 
gonal  to  the  contour  should  be  retained  from  in  the  shadow  in  the 

vicinity  of  the  contour  L.. 

In  the  future,  instead  of  the  function  U  (x,  x';  k)  it  will  be  more 
convenient  to  consider 


_  I 

I  Jf';  Uix,  A  ;  A). 


<2.S) 


By  distinguishing  the  reflected  and  enveloping  waves,  we  will  write 
A  (x,  x';  u)  and  (x,  x';  k). 

The  terms  with  second  derivatives  in  equation  (2,7)  should  be 
preserved  only  if  the  derivative^’  are  of  a  higher  order  than  the  function 

i 

U  (x,  x';  k)  itself  when  k  -►+  oo.  Then  V*  ^  is  of  principal  order 
in  the  expression  y  ^U;  therefore  instead  of  (2,7)  we  can  write  asymp¬ 
totically 

V*'A-{-iA|2v<I>vA-f  v-<l>Al=0.  ,,  o\ 


The  radiation  conditions  will  be  used  in  the  form 
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’  0,  /?  =- 1  I  ->■  . 

(2.10) 

Here,  (<}>)  is  the  arc  of  the  circle  |  xj  ■  R,  lying  in  the  region 
where  the  phase  <j>  is  defined.  By  using  the  asymptotics  of  the  Hankel 
function  when  <}>  -*  oo  and  the  Cauchy-Schwarz-Buniakowski  inequality, 
we  can  prove  that  the  radiation  condition  in  the  form  (1.  2)  follows  from 
condition  (2.10)  for  each  component  of  r((|>). 

2.  Let  us  examine  equation  (2.  9)  for  the  enveloping  wave. 

For  definiteness,  we  will  consider  that  a  ray,  directed  from  point 
x'  to  X,  envelops  contour  L  in  a  clockwise  direction.  We  will  also 
agree  to  measure  the  arc  length  clockwise. 

Let  us  use  two  orthogonal  coordinate  systems:  fs,  <)>g}  and 
{o',  n}  (figure  3). 


Figure  3  ^ 

Let  t  be  the  length  of  a  segment  of  the  tangent  of  Hx  and  s  the 
length  of  the  arc  at  the  point  of  tangency  H.  The  quantities  t'  and  s', 
associated  with  the  point  x',  are  introduced  in  a  symmetrical  manner, 
♦g  *  t'  +  (s  -  s')  +  t.  The  Lame‘s  coefficients  of  the  system{s,  <(>gjare 

~  TuT  ’  ^  * 


(2.11) 
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For  the  second  coordinate  system,  n  is  the  length  of  the  nor¬ 
mal  to  the  contour  L,  drawn  through  the  point  x,  a  is  the  arc  length, 
corresponding  to  the  point  where  the  normal  and  contour  intersect, 
n'  and  c'  are  defined  analogously.  The  Lame'^ coefficients  are 

=  '  +  -  *»=1.  (2- 12) 

Using  the  system  of  coordinates  [s,  ^g,}we  get 

=  (2.13) 


and  equation  (2.  9)  takes  the  form 


V I  =.0.  (2,  U) 

If  we  retain  here  only  the  term  in  brackets  containing  the  factor 
k,  we  will  arrive  at  the  equation 


with  the  boundary  condition  (2.4): 


•^0 

Solution  of  this  problem  is: 

A^^>(jc,  jc':  A):^0. 


(2.15) 


(2.16) 


This  geometric  optics  approximation  is  inadequate  for  our 
purposes.  To  construct  a  more  exact  asymptotic  of  the  enveloping 
wave,  let  us  examine  equation  (2.14),  at  first  in  the  vicinity  of  con- 
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tour  L  (the  else  of  this  neighborhood  will  be  defined  below).  Taking 
account  of  remarka  made  in  the  preceding  paragraph*  we  have 
asymptotically: 

Aa  +  «[2v0.,  vA''“  +  -;-  A''’]  =  0.  ^2.17) 

iS) 

Here  represent!  the  second  derivative  with  respect  to  the 

variable  n*  for  fixed  (r .  Similar  notations  are  used  later. 

Examining  equation  (2.17)  near  the  contour,  let  uo  use  the 
method  first  used  by  M.  A.  Leontovitch  and  V.  A.  Fok  [7]  in  diffrac* 
tion  problems  in  which  the  radius  of  curvature  was  constant.  This 
is  called  the  parabolic  equation  method.  V.  I.  Ivanov  [  3]  used  this 
method  for  contours  of  a  general  form.  We  will  refine  the  results 
of  the  work  (  3]  somewhat. 

In  the  equations,  the  differential  operations  and  coefficients 
will  be  represented  by  the  variables  n  and  <r  in  the  vicinity  of  the 
contour.  Let  us  use  the  formulas 

t~V2f(3)n’ - j-s'(a)n  +  o((i'‘). 


Let  us  retain  the  following  terms  in  equation  (2.17): 


(2.19) 


+  A'« + j.  A«.J  =0. 
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If  w«  conaidar  this  aquation  in  the  variahlee  ni 


n  and  r| 


■  r. 


it  ie  easy  to  see  that  terms  of  the  two  higher  orders  with  respect 
to  the  variable  k  (of  the  orders  and  k)  are  retained  in  equation 
(2,17)  and  the  discarded  terms  are  of  an  order  not  exceeding  k^  • 


Equation  (2.19)  should  be  considered,  therefore,  when 


1. 


Later,  when  analysing  (2.19),  we  will  not  discard  anything. 
Let  us  rewrite  (2. 19)  in  the  variables 


t  I 


t  I 


y  (4)*  J 


It  assumes  the  form 


■+4(4)'%a<-]=o. 


(2. 20) 


(2.21) 


(2.22) 


If  we  define  V'  '  by  means  of 


we  get 


(2.23) 


(2.24) 
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Finally,  let 


H.  K). 


(2.25) 


where 


^i(!‘ 


(2.26) 


then 


(2.27) 


In  definitions  (2.23)  and  (2.25),  the  symmetric  dependence  of 
(S) 

A'  '  (Xi  on  the  points  x  and  x*  is  explicitly  taken  into  account, 

where  the  point  x'  was  assumed  to  be  located  in  the  vicinity  of  the 

A 

contour  n'  ■  k*  n'  «  1.  Solutions  of  (2.27)  satisfying  the  condition 


,  0.  K)  =  0 


(2.28) 


and  which  have  a  derivative  quite  rapidly  approaching  zero 
when  p  «  (which  corresponds  to  the  radiation  condition,  since 
the  point  x  recedes  from  the  contour  along  the  normal  when  p  oo) 
can  be  represented  by  the  contour  integral  (when  p  ^  p') 


f',  =  f.  f  X 

X  W.  (C  -  lO  1  • 


(2. 29) 


Here  W|  (t)  and  (C)  are  Airy  functions  as  defined  by  V.  A.  Fok  [8). 
They  satisfy  the  equation  w"  (U  •  Cw  (C)  and  are  entire  functions  of 
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the  argurnciit  itscU.  '1  he  zeros  ol  the  I'unction  M](C)  are  located 
in  the  ray  arg  &  ‘  have  a  positive  imaginary  part.  Further 

properties  of  tliesi-  functions  can  be  found  in  [8],  The  integration 
contour  in  (2.29)  encloses  the  zeros  of  the  function  wj  (^);  the  fac- 
tor  e  ^  ensures  coiu  ergence  of  the  integral.  An  arbitrary  analy- 
tical  function  !•'( -j)  will  be  det«!rjnined  later  from  the  conditions  of 
merging  with  the  asyinj)tutic  in  the  light. 

The  formulas,  contaiJiing  the  integrals  (2.29),  were  obtained 
first  by  Fok  (see,  e,  g.  [2]),  in  the  case  when  p  (a)  *  const.  Later, 
they  were  examined  for  other  contours  admitting  separation  of 
variables,  e,  g. ,  for  the  ellipse  in  [4],  for  the  parabola  in  [  9]  And 
in  other  works. 

Formulas  (2,23)  -  (2,29)  give  the  asymptotic  in  the  vicinity 
of  contour  L,  It  can  be  assumed  that  during  a  suitable  change  of 
the  arguments  \  and  p,  this  fbtmula  will  also  give  the  asymptotic 
far  from  the  contour,  JLel  us  express  the  variables  \  and  p  in 
terms  of  t  and  s,  which  arc  more  natural  from  the  view'point  of 
the  geometry  of  the  problem.  Near  the  contour 

.'>(o)  P  t-s)  :  O  Cr  )  ; 

i  I=  ! 

'<0  .j,/.;.  (2.30) 

•  » 

■I  p  '  :  P  ■’  (V)r  }  Oin), 


In  (2.25),  let  us  sViballVute  the  new  variables 
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j*  r== 


fU^)  '  • 

V-'io*  • 


I  I 

i-  / 


<2.  31) 


In  (2.23),  let  us  substitute 


f.(a)^p(s). 


(2.32) 


In  the  vicinity  of  the  contour,  the  old  and  new  variables,  by  virtue 
of  (2.  30),  agree  to  an  accuracy  of  a  lower  order  of  magnitude  with 
respect  to  n(n  0). 

It  is  easy  to  prove  that  equation  (2.24)  for  the  function 
(X.  Y.  Y') 

. '1'  i  J’JK*  y  '  0  (2.33) 

can  be  written  as 


,  M*  (\> 

■'  fM.) 


(•') 


y^I».  -  -0. 


(2.34) 


(S) 

Here,  Vyy  represents  the  second  derivative  with  respect  to  Y  for 
fixed  X  and  Y',  and  the  gradient  is  calculated  with  respect  to  x  for 
fixed  x'.  Since  the  factor,  separated  out  by  relationship  (2.23)  (with 
the  change  <r  s),  is  carried  along  by  the  differential  operator 
V  into  equation  (2.  34),  then  (2.  34) means  that  in  the  initial 

X  d  X 


.1 


I  >  i  I 


(w.l  l)  till-  Si  t  Olid  c  uiiijiumjtil  ia  retained  coin- 


i  1  V  <  '  ‘  I 


•  ■  I  1  1  <•  1) 

.<■  t .»  ‘  “  ■ 

ih  arp.t  v'lit  I  d  vuil  oi  ;ln’  i id  lie  i,i  n  aper.ilor  V 

Wlieii  p  (a)  =  t,  oiiat,  tiu  i.  unattuu  y  <>l  X  me  iiib  that  (}i^  -  const;  in 
this  case  the  written  term  ia  to  the  second  derivative  of  the 

function  alon>;  tlu*  curvi.'  <{>_  =  conat.  The  curve  4)^  :  const 

intersects  the  geometric  light  -  sliadow  bouml.iry,  and  the  contour  L 
normally:  therefore  equation  {2.,  id)  makes  it  possible  to  assume 
that  in  the  variable  (2,  31),  the  asymptotic  formulas  W'ill  cancel  the 
principal  terms  o£  the  residual  (1.  3)  everyw'here  in  the  shadow.  When 
the  radius  of  curvature  is  variable  the  written  term  does  not  entirely 
correspond  to  the  derivatives  with  respect  to  the  direction  orthogonal 
to  the  boundary.  It  will  be  shown  in  section  3,  however,  that  the 
constructed  expression  nevertheless  cancels  the  principal  component 
of  the  residual  when  k  -♦  +  '«.  It  should  be  noted  that  in  the  vicinity  of 
contour  L,  expression  (2.23)  in  the  variables  \  and  p  more  accurate¬ 
ly  describes  the  asymptotic  since  it  cancels  the  two  higher  terms  in 
the  residuals  K  (x,  x';  k)  (1.3). 

3.  Let  us  turn  to  the  construction  of  an  asymptotic  expression 
for  (x,  X’:  k). 

An  orthogonal  system  of  coordinates  {  s,  (|>|^}  is  used  in  the 
light;  8  is  the  length  of  the  arc  at  the  point  of  reflection  P  (figure  2), 


IH  - 


"  t  +  t',  I  ia  tin;  h  nj;th  ul  ilu-  ai  j;im*nt  Px,  t'  Hu:  U  ngth  ul  «ignu‘nt 

x'P.  'I'hi'  I'orffic'irnia  ,»ri- 
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A 


(^.  ^b) 


0  thi:  angle  of  rofl  ti  tioj»,  /  ..  f  j(,\ 

It  is  i'asy  to  see  that 
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«4'A' 
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(2.  37) 


whore  on  the  right  side  of  this  <.'(|vtation(  i  is  understouti  as  a  func  tion 
of  8  and  Lot  us  write  the  geometric  optics  approximation 

^  is  usually  used  far  from  the  shadow  ,  It 

is  determined  from  the  equation 


-I- 


Oy) 

y 


€/'l 


0  (2.38) 


(compare  with  (2.9))  and  the  boundary  condition 

I  _>  j 

(CM.  2.1). 


thence 


we  got 


(‘2.39) 


(2.40) 


However,  using  (2.40)  in  the  light,  we  werVi  rtbl  able  lb  iielect 
the  function  P((,)  i*'  fhe  expression  fur  A  (see  (2.29)),  so  as  to 
obtain  condition  (2.  5)  of  continuity  of  the  asymptotic  at  the  boundary 
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KlV-T-<[vW  +  2lT^^^((  +  (.(s)cosO).|-iA<i  T=<t.;,V/<«lJ=:0.  (2.44) 

IR) 

The  subscripts  of  the  function  V'  ’  denote  derivatives  with  respect 
to  the  corresponding  arguments;  v**!*!^  given  by  formula  (2.37). 

The  coefficients  of  the  derivatives  in  relationship  (2.44)  can  be  ex> 
pressed  in  terms  of  the  variables  X,  Y  and  Y'. 

Let  us  put 


(2.45) 


where 


y.  y)=-^„x+-^(y-  :«)- + 


(2.46) 


is  determined  from  the  equation 


A'-(K-c,)”  -t-2(-;,)=  =0. 


(2.47) 


K  X,  Y  and  Y'  are  defined  by  the  formula  (2.41),  then 


(  —  Ca')'  =  A/  {s)  cos  0. 


(2.48) 


(R)  (R) 

4^*  '  is  considered  as  a  function  of  those  same  arguments  as  is  V' 

r  ^X,  Y,  Y',  s;  k).  Equation  (2.44)  takes  the  form 


(2.49) 
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The  last  relationship,  which  agrees  exactly  with  (2.27),  can  be  con- 

(R) 

sidered  as  the  equation  for  the  function  'I''  '  of  the  variables  X  and 
Y;  Y'  and  s  can  be  considered  fixed,  the  region  of  variation  of  the 
variables  X  and  Y  in  the  plane  {X,  Y)  is  included  between  the  two 
parabolas, 

[/T  +  APOO-  AI(s)]+  •C.Y< 

y/.’  (2.50 

(see  (2.47)).  It  is  necessary  that  the  function  (X,  Y,  Y',  s;  k) 

be  symmetric  with  respect  to  Y,  Y';  then  the  solution  of  (2.49)  can 

be  written  in  the  form  of  a  contour  integral  with  three  coefficient- 

functions  of  the  variable  of  integration,  of  the  variables  s  and  k. 

The  radiation  conditions  (2.10),  the  boundary  condition  (2.4)  on  the 
(R)  R  2  I 

contour:  ’Ir'  'e  X  I  y  _  q  =  -  1  requirement  that  the 

(R)  ”^^R  (R) 

function  V'  '  =  e  '  be  slowly  changing  (in  the  same  sense  as 

applied  to  this  term  in  paragraph  1  of  the  paper),  are  all  used  to 

( R) 

determine  these  functions.  The  function  '  of  such  a  type  (it  is 
written  when  Y'  ^  Y)  satisfies  both  the  radiation  and  boundary  con¬ 
ditions: 


I 
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O.) 


'' i‘- r,—  >") 


'  t .  I 


y)  . 

'•  I  Vs  J 


(2.  51) 
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The  second  integral  in  (2.  51)  with  the  arbitrary  func  tion  F  j(^)  vanishes 
on  the  contour  L,  The  first  integral  satisfies  boundary  condition  (2.4), 
which  is  easily  proven,  by  using  the  formula 


where 


.  V  -.i"'  ^  ( , 


(2.52) 


}■)-}- 


(2.53) 


V.  A,  Fok's  work  [10]  contains  a  proof  of  formula  (2.  52). 

Under  certain  general  assumptions  about  the  function  F  j  (4),  the 

second  integral  in  (2.  51)  can  be  computed  asymptotically  when 

according  to  the  saddle-point  method  (e.  g.  ,  [2]).  The  corresponding 

(R) 

component  in  the  function  V'  '  in  this  ease  will  have  a  rapidly  chang¬ 
ing  factor.  Therefore,  F|(4)  *  0  should  be  proposed. 

It  Is  seen  from  (2.  51)  that  depends  only  on  X,  Y  and  Y',  so 

that  later,  we  will  write  {X,  Y,  Y')  and  (X,  Y,  Y'). 

Let  us  determine  the  coefficient  F(^)  in  (2.29)  for  the  enveloping 
wave  so  that  the  continuity  condition  (2.  5)  of  the  asymptotic  in  the  geo¬ 
metric  light- shadow  boundary  will  be  fulfilled.  Beforehand,  let  us 

(R) 

transform  the  expression  for  V'  '  (X,  Y,  Y’)  by  using  (2.  52) 


Hero,  the  function  (2.29)  is  roprosonted  by  *  i 

F(t)  =  1.  The  difference  «(X,  Y,  Y')  -  Fj^(X,  Y,  Y')  can  be  re¬ 


presented  in  the  form 


(2.55) 


F(U  =  1. 


follows  from  formula  (2.  54)  and  the  merger  condition  (2,  6). 


(2.56) 


4.  Let  us  present  final  asymptotic  expressions  for  the  slowly 
changing  coefficients  (x,  x*,  k)  and  (x,  x';  k) 


where 


^  A';  A-)— A';i)  = 

a=A'''  V“'>(,Y.  V)  ^  y_ 


(2.57) 


(2.  58) 


The  function  F„  (X,  Y,  Y')  and  the  variables  X,  Y,  Y'  are 
Ja 

determined  by  formulas  (2. 46),  (2.47)  and  (2.41). 

For  the  enveloping  waves, 

i/'-'’  (X.  X';  k)  •-=  (.V',  A-;  k)  --  oj  (a,  a';  /c) 

('MjH'llr)'  J")  - 


where,  when  Y'  >  Y 
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(2.  60) 


The  contour  of  integration  encompaeseB  the  zeros  of  w  |  (C)«  The 
expressions  for  Fg  (Y#  Y')  and  the  arguments  X,  Y,  Y'  are  given  by 
formulas  (2.26)  and  (2.31). 

The  exponential  decrease  when  Z  -f  oo  is  an  important  property 
of  the  function  ^  “  Y^  +  Y'*  +  Z,  Y,  Y'  this  decrease  is  ob¬ 

tained  easily  by  drawing  the  contour  of  integration  such  that  Im  t  >0. 

If  the  errors  of  the  asymptotic  formulas  are  estimated  by  the  absolute 
value,  the  exponential  decrease  makes  it  possible  to  be  limited  to  a 
finite  number  of  enveloping  phases. 

Let  us  consider  the  asymptotic  Q(x,  x';  k)  which  includes  only 

/ 

such  phases.  When  the  points  x  and  x'  are  in  the  shadow,  there  are 
the  two  smallest  enveloping  phases  (figure  I).  When  x  and  x'  are  in 
the  light  relative  to  each  other,  the  phases  are:  direct,  reflected  and 
one  least  enveloping. 

The  asymptotic  Q  (x,  x*;  k)  has  properties  a,  b,  c  and  d,  entuner- 
ated  in  section  1,  Only  c  and  d  require  discussion.  Property  c 
follows  readily  from  the  uniform  continuity  of  the  functions  Y.  V) 

and  X*  (X,  Y,  Y')  in  the  octant  X,  Y,  Y'  >0.  The  functions 

(X.  Y,  Y')  due  to  the  exponential  convergence  of  the  integral  (2,  60), 
and  (X,  Y,  Y')  due  to  (2.54)  are  infinitely  differentiable  functions 


of  their  arguntents  when  X  >  0.  Thc-refore  with  fixud  x',  the  functions 
(R)  (S) 

U'*''  (x.  x';  k)  and  U'  '  (x,  x*;  k)  are  twice- continuously  dilf*.  rentiable 

with  respect  to  x  up  to  the  boundary  <j>g  •  This  means  that  the 

singularities  of  the  derivatives  of  the  function  Q(x,  x';  k)  agree  with 

those  listed  in  condition  d.  It  should  be  noted  here  that  the  component 

in  the  asymptotic  which  corresponds  to  tlie  smallest  enveloping  phase 

in  the  light  can  introduce  in  the  derivatives  discontinuities  of  the  first 

type  along  the  line  on  which  the  smallest  enveloping  phase  is  determined 

non-uniquely.  The  enveloping  wave  on  this  line  is  exponentially  small 

when  k -♦  4  oo;  therefore  it  can  be  corrected  such  that  the  derivatives 

become  continuous,  and  the  error  admitted  with  such  a  correction  will 

_k  3  A 

yield  a  contribution  to  the  residual  of  the  order  of  e  ,  where 
A  >  Aq  >0  and  Aq  is  not  a  function  of  x  and  x*.  Later  such  a  correc¬ 
tion  will  be  implied  without  explicit  mention  of  it. 


§  3.  AN  ESTIMATE  OF  THE  ORDER  OF  THE  RESIDUAL 


1,  The  residual  K  (x,  x';  k)  (1.3)  for  the  constructed  asymptotic 
Q  (x,  x';  k)  will  not  be  uniformly  small  when  k -*  +  oo.  For  equation 
(1.4),  however,  it  will  be  sufficient  to  restrict  ourselves  to  a  certain 
integral  estimate  of  the  residual,  for  example  to  the  following:  for 
any  continuous  function  p  (a)  (  a  >  0)  such  that  P  ^ 


n:ax  f  Uvpil: 

a,.:j 


-'••-vD./Vo'. -v;/;) 


/w  •  ►  'I*  CO  , 


(3.1) 
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The  complete  estimate  (3.1)  will  be  given  elsewhere;  here  we 
restrict  ourselves  to  an  estimate  of  the  order  of  the  residual  with 
respect  to  the  variable  k.  We  will  estimate  individually  the  expressions 


k)  and  —  *). 


and  also  the  contribution  to  the  residual,  which  comes  from  the  dis¬ 
continuity  of  the  derivative  of  the  asymptotic  Q  (x,  x';  k)  on  the  geo¬ 
metric  light- shadow  boundary.  These  expressions  can  be  estimated 
by  functions  such  as  M^(s)  f  (X,  Y.  Y')  (i  -  0,  1,  2)  (see  (2.  31)  and 
(2.41)).  Such  estimates  separate  the  dependence  on  the  variables  k 
and  X,  Y,  Y'.  With  respect  to  the  properties  of  the  functions  f  (X,  Y,  Y'), 
let  us  note  only  that  they  are  bounded  for  large  values  of  their  arguments. 
All  estimates  will  be  made  when  X  >  Xq  >  0;  this  restriction  is  not 
essential  for  integral  estimates  such  as  (3. 1). 

2.  Let  us  consider  the  contribution  to  the  residual  from  the 
iveloping  wave: 


( V ■  +  A’*)  r  («I. j  ic:  (v;;  4  k')  -J-  (.v,  a-';  k)  - 

x(2| K'M- 4 1"'']  i- i,;;  t;’  a‘^'+ 


(3.2) 


where 


p  <[ 
t  Os 


(3.3) 
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The  derivative  with  respect  to  t  will  be  taken  partially  relative  to 
a,  and  with  respect  to  s,  partially  relative  to  Formula  (i.2) 

is  obtained  from  the  last  expression  of  (2.  6)i  if  the  orthogonal  system 
'  coordinates  (s,  <>^3  is  used  and  if»  instead  of  the  partial  derivative 
\>  h  respect  to  4)^,,  the  derivative  with  respect  to  t  is  computed: 
for  £>  ■  const  ^  .  Let  us  note  that  k4>g  -  2M*  (8)x  ^  (s)  X, 

where  ^(s)  >  0  is  uniformly  bounded  from  above  and  below*  Writing 
the  differential  operations  in  (3.2)  in  terms  of  the  variables  X  and  Y, 
we  will  find  that  the  first  and  second  components  in  (3.2)  are  estimated 
by  an  expression  such  as  M  (s)  f  j  (X,  Y,  Y')  when  k  f  oo.  The  last 
component  seems  more  complex 


'  M  f  (a) 


(3.4) 


Here  Ls^^,  Lj^^  and  l\®^  are  universal  (independent  of  the  contour) 
operators  no  higher  than  the  second  order  with  respect  to  the  variables 
X,  Y  and  Y';  the  coefficients  of  the  operators  also  depend  only  on 
these  variables.  The  first  term  in  (3,  4),  of  a  higher  order  when  k  -•  +  oo, 
is  equal  to  zero  in  view  of  differential  equation  (2.  33).  Since  all 
terms  of  lower  order  are  proportional  to  the  derivative  of  the  radius 
of  curvature  p(s),  expression  (3.4)  vanishes  in  the  case  of  a  circle. 
Generally,  the  last  component  of  discrepancy  (3.2)  is  estimated  by  a 
function  such  as  M^(8)f2(X,  Y,  Y'). 
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Let  us  note  again  here  the  estimate 


+  A’) r (‘i'j  1  ‘Ac  • :  MX,  r;  r) dXdr, 


(3.5) 


which  uses  the  expression 


iix 


1 


dXdY. 


(3.6) 


for  the  volume  element.  The  right  side  in  (i.  5)  approaches  zero  as 

-  J. 

k  ^  when  k  -»+  »jo,  if  the  volume  element  dXdY  does  not  depend  on 

k. 

3.  The  residual  (r*  •'  A*)r(<!g  -(v;  -I-  A’):  ^  'j  //o ’(A^I^)  x';  k), 

which  it  is  also  convenient  to  write  as  the  sum  of  three  components 
(in  order  to  separate  k  and  X«  Y*  Y'): 

(ri+*>)i  ('iv)=  '-  (//;,"■  X 

+  /Al"(Mv)  4  {»;.  '■ 

is  estimated  analogously,  where 


(3.7) 


U  c _ *\ 

^  7  co.s  0 


(3.8) 


The  derivative  with  respect  to  t  is  taken  when  s  -  const  and  with 

respect  to  s  when  <)>^  ~  const.  When  deriving  (3.7)  from  (2.6),  the 

system  of  coordinates  (a,  <t>^]  is  used,  and  the  derivative  with  re- 

'  3  3 

spect  to  is  replaced  by  the  derivative  with  respect  to  t;  --gp 

R 
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(when  •  =  const).  The  relation  k«}»^  =  2M*  (s)  X  is  exact.  The  first 
two  components  of  (3.7)  have  an  estimate  M(8)  £4  (X,  Y,  Y')  (2.41). 


The  last  component  is  such: 

4  A-'  (.1  H-  i/'«  +  /•';;  ]  + 

I  . 

■  «->•+■'  -"i'r.Ti'A'''''-'!  /.fV'*'  1-  ,3.  „ 

*  ^  t  ..k  • 


+ 


(R)  (R)  (R) 

L3  ,  L2  '  and  L  j  '  are  the  differential  operators  of  an  order  no 
higher  than  the  second  with  respect  to  the  variables  X,  Y,  Y'.  In 
contrast  to  the  previous  case,  the  derivative  is  also  taken  with  re¬ 
spect  to  Y'.  The  higher  order  terms  are  cancelled  as  before.  For 
a  circle,  only  the  third  and  fourth  components  are  preserved  in  the 
braces,  and  the  last  component  of  (3.7)  is  of  the  same  order  as  the 
first  two.  For  a  general  contour,  the  estimate  is  M*  (s)  f^  (X,  Y,  Y'). 

Similarly  to  (3.  5) 


.  1 


(v;  +  r  (<!>,.)  j  <ix  •:  „  (.v.  r.  r)  c!X<ir 


(3.10) 


4.  The  discontinuities  of  the  first  order  derivative  of  the 
function  Q  (x.  x';  k)  in  the  geometric  light- shadow  boundary  reduce 
the  singularities  of  the  6-function  type  in  the  residual  K  (x.  x';  k). 
The  coefficients  of  these  singularities  are  proportional  to  the  jumps 
VjjQ  (x.  X*;  k).  In  the  operator  K,  applied  on  bounded,  continuous 
functions  by  the  expression 
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W){x)~  f  i1yf{y)K{y,A‘,  I:), 
h 


(3.11) 


the  last  component  corresponds  to  the  just  indicated  singularities 


iQ'n'in.  x-,k)-  Q''> -t; 


(3.12) 


The  integral  is  taken  along  the  geometric  light- shadow  boundary; 

(y^,  x;  k)  and  (y|i  x;  k)  are  the  limiting  values  of  the 

derivatives  at  the  boundary  in  a  direction  normal  to  it  from  the  shadow 
S  (x)  and  from  the  light  R  (x)  (positive  direction  of  the  normal  into 
S  (x)).  Instead  of  the  derivatives  with  respect  to  the  normal*  we  can 
compute  the  derivative  with  respect  to  the  lines  (J>g  =  const  and 
=  const. 

I'or  this  component,  (3.1)  denotes 


II 


r  if'  1  -  y, I  )I  Q;;''  O’,.  A':  k)  -  /,)|  ..q,  k  -h  co.  (3. 13 ) 


V  i.-iuicr  in  more  detail  the  jump  in  the  normal  derivative, 
:)  niiikes  it  possible  to  write  an  expression  for  the  sum 
>  t  .i!ui  reflected  waves  in  the  light  in  the  form 


I 


(3.14) 


The  variables  X,  Y  and  Y'  have  the  value  (^.41).  The  difference 
U  -  Fd  (see  55))  and  its  first  derivative*  with  respect  to  the  normal 

I\ 
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The  right  aide,  aa  before,  will  be  of  the  order  k  3  with 
dependent  of  k. 

The  author  thanks  O.  A.  Ladyzhenskaia  and  L.  D. 


dY  in- 


Faddeeva 


for  their  assistance  and  remarks. 
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